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■ Abstract We consider rotationally symmetric p-harmonic maps from the unit disk 

D 2 C M 2 to the unit sphere S 2 cl 3 , subject to Dirichlet boundary conditions and 
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with 1 < p < oo. We show that the associated energy functional admits a unique 
minimizer which is of class C°° in the interior and C 1 up to the boundary. We 
also show that there exist infinitely many global solutions to the associated Euler- 
Lagrange equation and we completely characterize them. 
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1 Introduction and main results 

Given a domain SI C K , a real number p > 1, a smoothly embedded compact 



submanifold without boundary M of ffi^" 1-1 , and a mapping u : Q i— > M, we consider 
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the energy functional 

Ep(u) ■=- [ |Vu| p dx. (1) 

The p-harmonic flow associated to E p is given by 

u t = -n u (-Apu), A p u = div(\Vu\ p ~ 2 Vu), (2) 

where, for a generic point v € M, we denote by n v the orthogonal projection of 
l^^ 1 onto the tangent space T V M to M at v. Notice that, because of this projec- 
tion, any solution to Eq. ([2J is constrained to remain in M. We next particularize 
M to be the unit sphere S N of R N+1 . Then, © can be written in an explicit form 

u t = div(|Vu| p_2 Vu) + u\Vu\ p . (3) 

This is a system of partial differential equations that, besides its pure mathematical 
interest as involving a competition between quasi-linear diffusion and reaction 
gradient terms, has been proposed in various contexts, such as ferromagnetism 
|10| , theory of liquid crystals [15] , multigrain problems [18] and image processing 
|22| . In this context it has been used as a prototype of quite complicated problems 
containing nonlinear reaction-diffusion systems of partial differential equations for 
the evolution of director fields. 

From the mathematical point of view, the Dirichlet problem for Eq. (|3]) with 
boundary condition u(t,x) = uq(x) for (t,x) £ dQ, Q := (0,oo) x ft, has been 
widely considered in the last decades. Referenced discussions can be found in [2] 
|3] for p = 2, [8 1TT1I20] for p ^ 2 and [13l[9l[TT] for the special limit case p = 1. 
About the local and global existence of solutions to the Dirichlet problem for Eq. 
([3]), it has been shown in [7] for p = 2 and, more recently, for p = 1 in |13|[TT] that 
for suitable boundary data, classical solutions exist for short time and their first 
derivative at the origin blows up in finite time. In [11] a sharp condition on the 
boundary data for this phenomenon to happen is given. In all these papers, the 
case of rotationally symmetric and stationary boundary conditions is considered. 

We devote the present paper to the study of the steady-states corresponding 
to the Dirichlet problem for Equation ([3]) when 1 < p < oo, M = S 2 , the unit 
sphere of R 3 , and Q = D 2 = {(x 1 ,x 2 ) £ K 2 : x\ + x\ < 1}, the unit disk of K 2 , and 
we restrict ourselves to deal with rotationally symmetric solutions. A rotationally 
symmetric solution of ([3]) has in this case the general form 

(X\ X"2 \ 2 

— sin h(t, r), — sin h(t, r), cos h(t, r) I , r = \x\,x = (x\, x 2 ) £ D (4) 
r r J 

and the Dirichlet boundary condition takes the form h(t, 1) = h(0, 1) = I. The 
energy functional E p introduced in ([1]) becomes 

E p (u) = / L p (r,h,h r ) dr, (5) 
Jo 

with a non-coercive Lagrangian given by 



L p (r,s,0 = (r 2 / p h 2 r +r ( 2 - 2 ^/ p sin 2 h) 
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We replace u by its special form ^ in Eq. ([3]). By a long but straightforward 
calculation, using the formulas given in the proof of Lemma 2.1 in [13] . we obtain 



A p u = (p~2)\Vu\ p ~ 4 h r 



h r h r 



sin ft h r sin ft cos ft 



(e cos ft, — sin ft) 



+ \Vu\v- 2 
'ft 



( -J + hr 



Sin ft + I h ftrr 

r 



cos ft 



+ ftr 



sin h — h r cos ft 



where we let 



(xi,X2) = re 16 . Thus, Eq. (J3j> becomes 
,2 h \ (p-4)/2 



ft, 



,2 . szn 



(p — l)h r h r 



+ CP-3) 



/i r sin /i cos ft ft r sin ft 



(6) 



— + 

r 



h rr sin ft sin ft cos ft 



Notice that, in the special cases p = 1 and p = 2, the previous equation simplifies 
to the equations in [9l p. 535], respectively [3l p. 96]. 

In the present paper we investigate the rotationally symmetric minimizers of 
E p subject to ft(l) = I, as well as its smooth critical points, that is, the steady 
states of Eq. ([6]). We thus fill in the gap between the limit cases p = 1 (studied in 
[9]) and p = 2 (studied in [3]), and moreover we extend our study to 2 < p < oo, 
giving at the same time an idea of how these results change with respect to p. By 
the periodicity of sin 2 ft and the symmetry with respect to 7t/2, we may assume 
without loss of generality that / £ (0, 7r/2]. We incorporate the boundary condition 
by defining 

(Ep(h), if = f, 

G p ,i(h) ■■= I (7) 
I +oo otherwise. 

We prove first existence and uniqueness of the minimizer of the energy functional. 



Theorem 1 Let p > 1 and I £ (0,7r/2]. Then, up to a reflection with respect to 
ft = 7t/2 in the case I = n/2, there exists a unique solution hi to 

ft z G argmin{G Pi; (ft) : ft G W hl (0, 1) CI W^(0, 1)} . (8) 

Furthermore, hi G C°°((0, 1]) n C 1 ([0, 1]), ft; is positive and increasing in (0,1], 
ft;(0) = 0, ft/(l) = I and hi solves the following ordinary differential equation 



( P l)ft-ft" + (p - 3) I ft' 2 Sinfc 2 C ° s/l - ft'^ 



h f3 h" sin 2 h sin 3 h cos h _ 



(9) 



in a classical sense. 
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Fig. 1 Global solution to (0 for 1 < p < 2 as given by Theorem [2] (a) and (b). Numerical 

simulation for p — 1 + "3 • Left: Global solution. Right: zoom of the global solution at the 
origin. 




Fig. 2 Global solution to (|9j for 1 < p < 2 as given by Theorem [2] (e) . Numerical simulation 
for p = 1 + Left: Global solution. Right: zoom of the global solution at the origin. 

We also show that the functional G p admits an infinite sequence of smooth 
critical points which are different from the minimizer. We distinguish between the 
cases 1 < p < 2 and p > 2. 

Theorem 2 Givenp £ (1,2), there exists a global solution h £ C°°((0, oo))nC 1 ([0, 00)) 
to Eq. (|9]) in a classical sense, which satisfies the following properties: 

(a) /i(0) = 0, lim h(r) = tt/2; 

r—>oo 

(b) There exists an increasing sequence of critical points r n of h such that r n — > 00 as 
n — > 00, h(r2n) £ (7r/2,7r) and are local maxima, h(r2n+i) & (0i 7r /2) and are local 
minima. Moreover, there exists some positive integer uq sufficiently large such that 
\h(r n ) — tt/2| is decreasing for n > no, and 

\h(ro) - || = max{|ft(r n ) - j| : n > o} ; (10) 

(c) For any I £ (0,7r/2), the minimizer h[ of Gp satisfies ftj(r) = h(ar) (or hi(r) = 
tt — h(ar) if I = tt/2), for some suitable a > 0; 

(d) Any non- constant finite global solution h to Eq. Q has the form h(r) = kn±h{ar) , 
for some a > and k £ Z. 

(e) There are infinitely many global solutions which have —00 (respectively +00) as 
limit when r — > + . They are either increasing (resp. decreasing) forever and with limit 
either +00 (resp. —00) or kir for some k £ Z as r — > +00 or they increase (resp. 
decrease) up to some r* £ (0,+oo) and after r* they oscillate around (2k + 1)-| for 
some k £ 1 as in (b). 

In order to illustrate graphically these solutions and ease the understanding of the 
results, we plot them in the following figures. 
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Due to the scaling invariance of Eq. (|9j with respect to r, Theorem [2] shows 
that, for 1 < p < 2, in addition to the minimizer given in Theorem [T] and to the 
constant solution h = tv/2, there exists an infinite sequence of critical points for 
G p subject to h(l) = tt/2. This is similar to the limit case p = 1, as proved in [9], 
but in strong contrast with harmonic maps; i.e. p = 2, where the minimizer is the 
only global monotone solution connecting to tt/2. 

For the case p > 2, the picture is completely different and it is surprisingly 
very different from the case p = 2. 

Theorem 3 Letp>2. Then there exists a global solution h £ C 00 '((0, oo^nC 1 ([0, oo)) 
to Eq. © in a classical sense, which satisfies the following properties: 

(a) h(0) = 0, h is increasing and lim h(r) = +oo; 

r— >oo 

(b) For any I £ (0,7r/2), the minimizer hi of G p satisfies /i;(r) = h(ar) (or hi(r) = 
tt — h(ar) if I = tt/2), for some suitable a > 0; 

(c) Any non-constant global solution h to Eq. Q such that h(0) = kn for some k G Z, 
has the form h(r) = kn ± h(ctr), for some a > 0. 

(d) There are infinitely many global solutions to ([9]). They verify h(0) = ( 2fe + 1 ) 7r j or 
some k £ Z. They may oscillate around ( 2fc + 1 ^ U p i some r* G (0,+oo) and then 
either they increase up to 2(k + f )7r or to +oo or they decrease up to 2kix or to — oo. 
fe) There are infinitely many global solutions which have — oo (respectively +oo) as 
limit when r — > . 27ie?/ are either increasing (resp. decreasing) forever and with limit 
+oo (resp. —oo). 

We plot the most interesting solutions given by Theorem [3] in the following figures. 



400 




Fig. 3 Global solution to J9]l for p > 2 as given by Theorem [3] Numerical simulation for 
p = 3. Left: global solution corresponding to case (a). Right: global solution corresponding to 
case (e). 




Fig. 4 Global solution to J9j for p > 2 as given by Theorem \3\ (d) . Numerical simulation for 
p = 2.1. Left: global solution. Right: zoom of the global solution at the origin. 



6 



Razvan Gabriel Iagar, Salvador Moll 



The paper is organized as follows: first, in Section [2] we show the existence 
(and regularity) of a minimizer through variational techniques. Then, we restrict 
to solutions of Eq. © and classify them, in order to prove our main results. This 
is done by the analysis of a rather complicated phase plane associated to ([9|) after 
various changes of variable. In a first step, we study in Section [3] the behavior of 
solutions h of ([9]) as r — > 0, which allows us to complete then the proof of Theorem 
[1] by showing the uniqueness of monotone solutions with the expected behavior 
for ([9]); the latter is done in Section [U In a second and final step, in Section [5l we 
prove Theorems [2] and [3] by analyzing the behavior of solutions of Eq. ([9]) for large 
r. 



2 Existence of a smooth minimizer 



First of all, we collect some properties of solutions to problem ([5]). We state the 
next result without proof since its proof is identical to that of [9] Lemmas 1—3]. 



Lemma 1 Up to a reflection with respect to hi = ^ if I 
problem 0) is nondecreasing and /i;([0,l]) C [0,1]. 



any solution h[ to 



Next, we show that any solution of (J5J) must be positive in (0, 1]. 

Lemma 2 Let h t G W M (0, 1) n W^((0,l]) with fy(l) = I be a solution to ©. 
Then, hi is positive in (0, 1]. 

Proof By Lemma [T] we assume without loss of generality that hi e [0,1] is nonde- 
creasing. Assume by contradiction that < a < 1 exists such that ft; (a) = and 
hi(r) > for a < r < 1. Let then e > such that a — e > and a + e < 1. We first 
notice that 



(i 



((hi) 



2+ sin^ 



dr > 



(■a-\-e 



rUhArfdr > min f r(h r Y dr, (11) 
heW J a 



with 
W 



| ft G W 1,p (a, a + e) : h nondecreasing, h(a) = 0, h(a + e) = hi (a + s) i 



We next estimate the last integral in (|11[) . In this case, the corresponding La- 
grangian is coercive and convex. Therefore, (see e.g [SI Theorem 3.7]) there exists 
a minimizer h a G W l,p (a,a + e). Moreover, it attains the boundary values and it 
verifies the corresponding Euler-Lagrange equation: 



r((ha)r) 



p-1 



C 



h a (r) 



1 2 — p 

(1 - p)C"- 1 r 1 -p 



+ D. 



Substituting with the boundary values we get: 



C 



{p-2){h[{a + E))[{a + E)p-i -qp-i] 



p-i 
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Then, 



mm 

hew 



a+e i>a-\-£ 
r(h r ) p dr= / 

J a 

(h l (a + s)n P ~2)P- 1 
(p-l)f- 1 



p -i 
C p- 1 rp- 1 dr 



(p- l)Cp- 1 rv- 1 



a+e 



p — 2 p — 2 _ p — 2 p — 2 

((a + e)p~ 1 - ap- 1 ) p ((a + e) p- 1 - a?- 1 ) 



p-1 



On the other hand, if we take 



9{r) = 2 arctan 



t an (M£±£l) 
(a + e) 



we have 



ra-t-e > 

/ Jr 2 (6") 2 + sin 2 6>dr 
Jo 

2f (a + £ ) 2 -P tan (M£±£)) P J cog 2p ( M^£) ) _ cos 2 



(1-p) sin 2 (M^£l) 



Therefore, letting 



we obtain that 



69(r), <r <a + e, 
hi(r), elsewhere, 



E p ,i{h*) - S PiJ (ft,) < ftf (a + e)(a + e) 2 " p x 

ip-i 



2-p 

(2 - pjop- 1 



.(p- l)((a + e)p-i -a"- 1 ) 
as £ — > 0. Finally note that 

(2 - p)^- 1 



+ 2~» +o e (l) 



(p- l)((a + e)»-* -op-') 



p-l 



+00 



for p ^= 2. Hence ft/ cannot be a minimizer. 

We next show that any solution to problem ([§} is smooth in the bulk. 

Lemma 3 Let hi € W 1 ' 1 ^, 1) n W^((0,1]) iratt ft,(l) = I be a solution to ©. 
Tften, ft, 6 C*°°((0,1]). 
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Proof In order to improve the regularity of hi, we claim that for any e > 0, hi is a 
minimizer of 

nl 

L p {r,h,h r ) (12) 



i: 



within the space of functions h such that h G W 1,p (e, 1) , h{e) = hi(e),h(l) = 
I, h nondecreasing. Let 

Kht : = { h G Wl ' P ( e ' !) 'Me) = = ^ • 

Then, since the Lagrangian in (JT2J is superlinear, coercive, analytic on (e, 1) x 
[hi(e), +oo) x R, convex with respect to /i r and with bounded derivatives with 
respect to h and h r (by Lemmas [1] and [2]) , then (0 Theorems 3.7 and 4.1]) there 
exists h e G W £ * such that 

ft e G argmin | ^ ft, for) : ft G ( e > ^) ' ^ nondecreasing | . 

Moreover, any minimizer is of class C°°([e, 1]). Observe that, once the claim is 
proved, then we obtain that 

ft; £ iv 1,1 ^, 1) n C°°((0, 1]). 

Suppose by contradiction that the claim does not hold. Then, defining 

' hi(r), if0<r<e, 

h e (r), if e < r < 1, 



ht(r) 



we would obtain that h* G W M (0, 1) n W^((0, 1]) and 



E p (h%) = L p (r,h*,(h* r ))dr + J^ L p (r, h £ , (h e ) r ) dr < E p {h{) , 

which contradicts hi verifying (|8|). 

These previous Lemmas permit us to prove the following existence result. 

Proposition 1 There exists a solution h of © such that h G W rl ' 1 (0, l)nC°°((0, 1]). 

Proof Since the Lagrangian is not coercive in [0, 1], the existence of the minimizer 
does not follow from classical results in the literature. We follow the direct method 
in the calculus of variations. 

Let {h n } C W ' p (0, 1) be a minimizing sequence such that /in(l) = I- By 
Lemma [TJ we may assume that each h n is nondecreasing and h n {[0, 1]) C [0,1]. 
Then, 

|| || M^ 1 . 1 (0,1) = ll^nlll + ll(' l «)''l|l 5= 2/. 

Thus, we can extract a subsequence, not relabeled, such that h n — > h* G BV(0, 1) 
in i 1 (0, 1). We also obtain that h n — > h a.e. in [0,1]. Moreover, given an e > 0, 
since _ 

"(A*)?dr<- / L(r,A„,(An) r )dr<^^<+oo ) 
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then there is a subsequence of the previous subsequence such that h 
W 1,p (e, 1) and /i**(l) = I. Therefore, since the whole sequence converges to h* , we 
obtain that h* € BV(0, 1) n W^((0, 1]) and h*(l) = I. Note that this implies that 
h* £ W 1>:1 (0,1) n Wf£ ((0, 1]). Moreover, since the Lagrangian is a Caratheodory 
function, then 

/ L(r, h* , h*) dr < lim inf / Z/(r, /i n , (h n )rdr < liminf E p (hn)- 
y £ n—¥oo J £ n— »oo 

Finally, by the monotone convergence theorem we conclude that h* is a minimizer 
since 

E p (h*) < liminf E p (h n ). 

n—too 

We finish the proof by applying Lemma [3] 



3 Local behavior at r = 

We begin with the systematic study of the steady states given as solutions to 

In the present section we deal with their behavior near the origin, in both cases 

1 < p < 2 and p > 2. 



3.1 The case 1 < p < 2 

The main result in this Section is the following one. 

Proposition 2 Let p £ (1,2) and h be a nonconstant solution to ([9]) in an interval 
(0, e) for some e > 0, and such that h(r) £ (0, n) for all r £ (0, e). Then h £ C 1 ([0, e)), 
h is monotone in a right-neighborhood of r = and 

lim h(r) = or lim h(r) = ir. (13) 
r->0+ r->0+ 

The proof of this proposition will be the consequence of some suitable changes of 
variable and a phase plane analysis of the resulting algebraic autonomous system. 
The changes of variable coincide with the ones in the particular case p = 1 [S]. 
First, we pass to logarithmic coordinates by letting f(t) := /i(e - *). Then, ([9]) 
becomes 

(P - l)/' 2 (/" + /') + (P - 3)[/' 2 sin /cos / + /' sin 2 /] - f' 3 

9 S ( 14 ) 

+ (/" + /') sin 2 / -sin 3 /cos/ = 0. 
For / £ (0,7r), as supposed in Proposition [2] we make the further change 

w(t) := cot f(t), 

in order to eliminate the trigonometric terms. We obtain the following algebraic 
differential equation 

„ _ (1 + w 2 + w' 2 ){2(p~ l)ww' + (2 - p)(l + w 2 ))w' 
W (l + w 2 )(l + w 2 + (p-l)w' 2 ) W ' { ' 
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which can be written in an equivalent form as the autonomous system 
w = k, 

, = (1 + w 2 + fc 2 )(2(p - l)kw + (2 - p)(l + w 2 ))k _ ^ 
{l + w 2 ){l+w 2 + (p- l)k 2 ) ™' 



The system (|16|) has a unique critical point in the plane, the origin. The linearized 
system around the origin has the matrix 



A(0,0) 



1 

-1 2-p 



Thus, the point is an unstable node, the matrix yl(0,0) having two complex eigen- 
values with real part (2 — p)/2 > 0, as it is easy to check. Then, all the orbits 
of (|16p come out from the origin and will connect to critical points at infinity. 
In order to study the critical points of (|16p at infinity, we note that (|16|) is not 
polynomial, thus the standard theory cannot be used directly. Instead, we write 

dk_ _ (1 + w 2 + k 2 ) [2(p - l)wk + (2 - p)(l + w 2 )] k 
dw~ k(l + w 2 )(l + w 2 + (p-l)fc 2 ) 

_ w(l + w 2 )(l + w 2 + (p - l)fc 2 ) _ Q(w,k) 
~ fe(l + w 2 )(l + w 2 + (p - l)fc 2 ) ~ P(w,k)' 

where Q(w,k), P(w,k) are both polynomials of degree 5. Thus, we can follow the 
strategy indicated in [21] pp. 270-271]. We let 

P(w,k) = Pi(w,k) + . . . + P 5 (w,k), Q(w,k) = Q 1 (w,k) + ... + Q 5 (w,k), 

with Pj, Qj homogeneous polynomials in (w,k) of degree j, 1 < j < 5. Then, 
Theorem 1 in pit p. 271] ensures that the critical points at infinity occur at the 
polar angles 9 solving the equation 

cos 6Q 5 (cos 6, sine) - sin 9P 5 (cos 6, sin 9) = 0. (18) 

Let us notice that, in our case, 

Q 5 (w, k) = k{w 2 + k 2 ) \2{p - l)wk + (2 - p)w 2 j - w 3 (w 2 + (p - l)k 2 ), 
Ps(w, k) = kw 2 (w 2 + (p — l)k 2 ). 
Hence (1181) becomes 



(19) 



cos # sin # |^2 (p — 1) sin 9 cos S + (2 — p) cos 2 

- cos 4 9 [cos 2 9 + (p - 1) sin 2 <?] = cos 2 sin 2 [cos 2 + (p - 1) sin 2 

We thus get six critical points at infinity, corresponding to polar angles 

. I 7T Sit tt 57r / p — 1 \ / p — 1 \ 

^ £ \ tt . "TT > T i > arccos = , - arccos = > . (20) 
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Since w = k, we cannot have at the same time k — > oo and w — > — oo or w — > oo 
and ^ -oo as i -> oo, that is, the product sin cos must be nonnegative at 
infinity. We thus reduce, also by symmetry, our analysis to the study of the two 
points corresponding to polar angles 9 = ty/A and 9 = tt/2. 

We start with the local analysis of the system (|16p near the critical point with 
9 = 7r/2, which is more involved. By Theorem 2 in [211 pp. 272-273], the flow in a 
neighborhood of this point is topologically equivalent to the flow near the origin 
(a;, z) = (0, 0) of the following system 

±a/ = zz B Q(zl)- x *p(*,l 
\z z j \ z z 

, / &r, fx 1 
± z = z Q — , - 

V z z , 

the signs plus or minus being determined by the flow on the equator of S* 2 as 
indicated in Theorem 1, [211 Section 3.10]. After straightforward calculations, the 
system above can be written in the form 

± x = (v — l)(x 2 — z 2 ) + higher order terms, 

t _ (21) 

± z = 2(p — \)xz + higher order terms. 



We will use next the theory developed in |19j to show that indeed the system ([21 
is topologically equivalent locally near (x,z) = (0,0) to the homogeneous system 



±x =(p - l)(x 2 - z 2 ), 
± z' = 2(p — \)xz. 



(22) 



In order to do this, we begin by passing in (|22|) to polar coordinates (r, 9) so that, 
after easy calculations, the system becomes 

dr 2,/^ d9 . . 

where 

f(9) := (p-l)cos0, 5 (6?) := (p- l)sin0. 

Introducing a new time s via ds/dt = r and making a further change of variable 
g := r/(l + r), we transform (|22[) into the topologically equivalent system 

±g = (p- 1M1- g)cos0, 
±6»' = (p- l)sin6», 

where derivatives are taken with respect to s, and 

(g,0) £fl:= {(g,9) : < 0< 1,0 G [0,2tt)}. 

This system has two critical points on dD, namely (g, 0) G {(1,0), (l,7r)}. The 
linearization of (|23|) near these two critical points has the matrices 

"<'•»>= K" " „-!))■ "('•') -Co' -0.-1)) ' 
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thus both points are hyperbolic (two saddles). We are then in the situation de- 
scribed in Theorems A and C of [19] , and from the latter theorem we deduce that 
(|2ip is topologically equivalent near the origin to (I22p , which is furthermore topo- 
logically equivalent, through a simple time rescaling r = (p— l)s, to the following 

, / 2 2 

± X = X — Z , 

+ ' 9 (24) 
± z = 2xz. 

It is easily seen that the origin is an elliptic point (see [1] for a proof). Coming 
back to the original system, it follows that the orbits coming out or entering the 
point with 8 = tv/2 describe elliptic cycles. 

We study now the second critical point at infinity, with polar angle tt/4. Since 
this point is hyperbolic (as we shall see), we pass to the Poincare sphere using an 
equivalent approach used for example in [TB] . We let 

w = — - — cos<P, k = — - — sin<P, 



whence 



Vw 2 + k 2 ., 1 _ k 

1 , , <P = arctan ■ 



1 + Vw 2 + k 2 1 + V™ 2 + k 2 w 

By a direct differentiation, we find that 

g ' = (1- qY(w cos <P + k sm<P) , $' = — = p3-. (25) 

w + £r 

Then, we arrive to the following equations 

' , _ e(l-g)[(l-e) 2 +e 2 ][2(p-l)e 2 cos<f sin*+(2-p)((l-e) 2 + e 2 cos 2 ^)]sin 2 * 

Q ~ l(i-e) 2 +e 2 cos 2 <2>J[(i-e) 2 +e 2 cos 2 s>+( P -i) e 2 sin 2 s>j ' 

< 

, _ [(l-e) 2 +e 2 ][2(p-l)£> 2 cos*sin#+(2-p)((l-e) 2 + g 2 cos 2 *)] sin^cos* 

, * — [(i-e) 2 +e 2 cos 2 <i>J[(i-e) 2 +£ 2 cos 2 &+(p-i)e J sin 1 4>\ 



(26) 



In this setting, our critical point becomes Pi = (l,7r/4). The linearization of the 
system (|26|) near Pi has the matrix 

-('■i)-(-o'i 

thus this point is a saddle. There exists a special orbit of (|16H going into Pi and 
coming from the interior of the plane. All other orbits are included in the boundary 
g = 1 and are not of interest to us. The orbit entering Pi and coming from the 
plane do this in infinite time and its behavior is 

lim HiM = 1 U m W M = oo. (27) 

t->oo w(t) t->oo 

The solution corresponding to this orbit is increasing for t > to large and it satisfies 

m 
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Proof (Proposition 0) We analyze the special orbit entering P± in the previous 
phase plane. Since lim w(t) = oo, then, exchanging if necessary h by iv — h, we 

t— f oo 

may assume that lim h(r) = 0. To complete the proof, it suffices to show that 

r— >0 

there exists lim h'(r) G R. 

In the following step, we want to prove that, along this orbit, we have w'(t) = 
k(t) > w(t) for t sufficiently large. To do this, assume by contradiction that there 
exists some e > and some to € (0, oo) such that w(to) = (1 + e)k(to)- Define 

l(t) := w(t) - {l + e)k(t). 

Thus l(to) = 0. Using (|16|) and the equality w(to) = (l+e)k(to), by straightforward 
calculations we obtain (everything being taken at t = to that we omit from the 
notation for simplicity) 



Z' = fc[(l + £ ) 2 + l-(l + e) 

[(1 + k 2 + (1 + e) 2 k 2 )} [2(p -!)(! + e)k 2 + (2 -p)(l + (1 + e) 2 k 2 )} 



[(1 + (1 + e) 2 fc 2 )] [1 + (1 + e) 2 fc 2 + (p - l)fc 2 ] 

W(k,p,e). 



(28) 



In order to bound from below \P(k,p, e) above, we first notice easily that J^&^kjp, e) > 
0, whence 

p, e) > 9{k, l,e) = l+e + e 2 - 1 ^^ 2j(2 - 

But the right hand side of the last inequality is decreasing with k, therefore we 
can take the limit as k — > oo and get 

e) > 1 + e + e 2 - jL- = > 0, 

which implies 

i'(t ) > ^M|+|£+£!) > o, 

hence w(t) > (1 + e)fc(i) for any t > to- This is a contradiction with the fact that 
w(t)/k(t) -» 1 as i -> oo. Thus fc(t) = u/(t) > tu(t). 
We next follow as in [9] by writing 

, // n rh'(r) tanft(r) 2,/ \ /™n 

h'(r) = . , f , y ' u , y -L C os 2 h{r). 29 

sin h(r) cos h(r) r 

On one hand 

rh'M ,• /'(*) ,• w'(t) , 

hm . ,,,,,, = lim — : — r ,V ir? x = lim — hr = 1. (30) 

r _>.o+ sm/i(r) cos/i(r) t->oc sm/(i) cos f(t) t-yoo w{t) 

On the other hand, we have 

^- t (w(t)e- t ) = e- t iw , (t)-w(t))>0, 
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for any t sufficiently large, by the previous step. This implies 




■Md = lim 7 \ _. = Ce [0,oo). 
r t— yco wltje 



(31) 



From pi?]). (|3"0|) . (|5T]| and the fact that ft(r) -> as r -> 0+, we deduce that 
ft'(r) — > C > as r — > + , which shows that h £ C 1 ([0,e)) for some e > small. 
The monotonicity in a right-neighborhood of is now obvious from the previous 
analysis. 

3.2 The case p > 2 

In this case, the behavior near the origin is slightly more complicate. More pre- 
cisely, we prove 

Proposition 3 Let p > 2 and h be a nonconstant solution to §§§ in an interval (0, e) 
for some e > 0, and such that h(r) £ (0,7r) for all r £ (0,e). Then h £ C 1 ([0,e)), h 
is monotone in a right-neighborhood of r = and 



Proof We do the same changes of variable as in Subsection 13.11 until arriving to 
the system (|16[) . We study the phase plane associated to (|16[) : the critical points 
are the same as for 1 < p < 2, that is, a unique critical point in the plane (the 
origin) and the six critical points at infinity corresponding to the polar angles 6 
given in (f20)l . 

Since now 2 — p < 0, the origin becomes in this case a stable node. By standard 
theory, this point is asymptotically stable [211 Section 2.9], that is, there exists 
a ball B(0,S) for some S > such that all the orbits entering B(0,S) end up at 
(w, k) = (0, 0). The local analysis of the critical points at infinity is the same as in 
Subsection l3.ll By performing the same analysis as in the proof of Proposition [2l 
we find the following two different types of orbits: 

• either the orbit entering the saddle point Pi of polar angle 9 = n/A and coming 
from the plane, that in initial variables means lim h(r) £ {0, 7r} as in the proof 



of Proposition [2] 

• or the orbits entering the attractor O, meaning that w(t) — > as t — > oo . In initial 
variables, it means that lim h(r) = n/2, ending the proof. 



4 Uniqueness of the minimizer 

We prove in this section the uniqueness of the smooth minimizer. More precisely: 

Proposition 4 Let I £ (0, 7r/2]. Then there exists a unique positive, non-constant and 
non- decreasing solution to the problem 



lim h(r) = or lim h(r) = — or lim h(r) 



(32) 



r->0+ 



r->0+ 



(p- l)h l2 h" + (p-3) K 



sin h sin 3 cos h 



,2 sin h cos h 




{ 




= 0, re (0,1), 



(33) 



r 



, K 1 ) = L 
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Furthermore, h is increasing. 

Proof The proof will be divided into several steps. 

Step 1. Monotonicity. We pass again to logarithmic coordinates as in the previous 
section by letting /(t) = h{e~ l ). Then, (|33f) transforms into (|14p posed for t € 
(0,oo), with boundary conditions 

/(0) = I 6 (0,7r/2], lim/(t) = 0, 

t— >oo 

with / positive, non-constant and non-increasing. We next show that, under these 
conditions, any solution / to (|14|) is in fact decreasing. Assume by contradiction 
that there exists to € (0, 00) such that /'(to) = 0. If /(to) = tt/2, by standard ODE 
arguments of uniqueness, we have / = tt/2, contradiction. Thus < /(to) < ir/2. 
Replacing t = to in (fTl)l . we get 

/"(to) = sin /(t ) cos /(to) = isin2/(t ) > 0, 

therefore / is increasing in a small right-neighborhood of to, contradiction. Thus 
f (t) < for any t £ (0, 00), and coming back to initial variables, it follows that h 
is increasing. 

Step 2. New change of variable. For / G (0,7r/2), we perform again the change 
of variable 

w(t) :=cot/(t), 

thus arriving to Eq. (|15[) . with boundary condition w(0) = cot/. Since we know 
from the analysis in the previous section that w is monotone for t sufficiently large, 
we can define its inverse t(w) such that w(t(w)) = w, and we make the further 
change of variable 

g ( w ) ■= w (34) 
yv ; w'(t(w)) v ' 

Differentiating with respect to w, we get 

„ _ g(w) - wg'(w) _ wg(w) - w 2 g'(w) 

^ W>> - g( w )H'(w) ~ g{w)3 

We substitute these formulas in Eq. (|15p and obtain 

g(w) - wg'(w) 

[g(w) 2 (l + w 2 ) + w 2 } [2(p - l)w 2 + (2 - p) g (iu)(l + w 2 )} g(w) 3 (35) 

(1 + ™ 2 ) [<?H 2 (i + w 2 ) + ( P - \)w 2 } 9[w) ' 

After performing some technical operations, we write (|35[) in a more suitable form: 
g'M = ^ [gH 3 - (2 - p)g(w) 2 + g(w) - F(w, g(w)j\ , (36) 



where 



_ w 2 g [2(p -!)(! + w 2 )g 2 + (2 - p) 2 (l + w 2 )g + 2(p - 1) W 2 ] 

(l + w 2 )[ ff 2 (l + w; 2 ) + (p-l)w; 2 ] ' 1 j 
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Step 3. Comparison estimates. Assume by contradiction that uniqueness as 
stated in Proposition U does not hold and pick two different positive solutions g\ 
and g 2 of (|36p . Adapting ideas from the case p = 1 [H] , our next goal is to estimate 
the difference g\ — g%. More precisely: 

(91 -92)'(w) = — \gi (w) 3 - g 2 (w f - (2 - p)( gi (w) 2 - g 2 (w) 2 ) 
{gi -S2)0) ■ 
(ffi -92){w) 



+ {gi - 92){w) - F(w,g 1 (w)) + F(w,g 2 (w))] 

2 



+(si32)(w) + 1 



9lW - 
(2-P) 2 



2-p 



+ 92(w) - 



2-p 



[F(w,gi(w)) - F(w,g 2 {w))] . 



Taking into account that lim gi(w) = lim g 2 (w) = 1, by the analysis in Section 

w— >oo w— >oo 

|3] we further obtain 



(gi - gi) 1 M 
{gi - 92)(w) 

' gi(w) 



2-p 



+ 92(10) 



2-p 



+ (g ig2 )( W ) + l-^y 



(38) 



_ 1 F(w,g 1 (w)) - F(w,g 2 {w)) > C _ 1 F(w,gi(w)) - F(w,g 2 (w)) 
w (gi-g2)(w) ~w w g 1 (w) - g 2 (w) 

for w sufficiently large, where the constant C satisfies 



2p. 



(39) 



In order to proceed with the integration of the differential inequality (|38|) , we only 
need to estimate conveniently 

F{w,g 1 (w)) - F(w,g 2 (w)) 

gi( w ) - g2(w) 

This will be our last step in the proof. 

Step 4. Lipschitz estimate and end of proof. In order to estimate the last 
quotient, we differentiate F defined in (|37p with respect to its second variable g: 

dF(w,g) _ 6(p - l)w 2 (l + w 2 )g 2 + 2(2 - p) 2 w 2 (l + w 2 )g + 2(p - l)w 4 g 
dg ~ (1 + w 2 ) [(1 + w 2 )g 2 + (p - l)w 2 } 

2g [2(p - l)w 2 (l + w 2 )g 3 + (2 - p) 2 w 2 (l + w 2 )g 2 + 2(p - l)w 4 g] 
[(l+«; 2 ) 5 2 + (p-l)«; 2 ] 2 _ ' 

Letting w — > 00 and recalling that lim g(w) = 1, we can compute 

w— >oo 

dF(w,g{w)) 



lim 

w— f 00 



dg 



2(P-1). 



(40) 
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Fix e € (0, 1). Then, by (|40[) we can find wq = wo(e) large enough such that 



-M 

31 ~92 

for any w > iuo(e), and at the same time 



<2(p-l)+e, 



cm 



3i(w) 



- p 



+ 



32 (w) 



2-p 



+ (3i32)0) + 1 



(2-P) 5 



> 2p- 



also for any w > u>o(e), the latter resulting from (|39[) . Integrating now (|38|) on 

(wo(e),iy), we obtain 



3i (w) - 92 (w) 





c 


w 






exp 


Wo 





. r 1 F (ai 

Jwn W 9 



) - H92) 



31 (wo) - 32 (wo) L w oJ * L Jw w 91-92 
for any w > wo(e), where 

C := 2p - 2(p - 1) - 2e = 2(1 - e) 
Since e < 1 by its choice, we get C > 0, whence 





w 


> 









(41) 



lim 



31 ( w ) ~ 92{w) 



00, 



w->-cx) g\(wo) - 32(^0) 

which is a contradiction with the fact that g{w) — > 1 as w — > 00, for any p solution 
of (|36p . This completes the proof. 

Proof (Theorem]!^ By Propositions[T]and[2]for 1 < p < 2, respectively Propositions 
[I] and [3] for p > 2, and Lemmas [T] and [2] there exists a solution /i to problem (P), 
and furthermore h € C°°((0, l])nC 1 ([0, 1]) with h(0) = and h(l) =£,his positive 
and nondecreasing in (0,1], and h solves Q. In order to show that h is unique, 
let h be any solution to problem ©. By Lemmas Q] [2] and |3] h G C°°((0, 1]), 
is positive and non-decreasing in (0,1], and h is non-constant since h(0) = 0. 
Hence, by Proposition |H h = h. As a consequence ft, is increasing, and the proof is 
complete. 



5 Smooth critical points 

5.1 The case 1 < p < 2 

In this subsection we prove Theorem [2] Proceeding along the ideas of [9], we study 
the behavior of the solutions h to Eq. © for r large. To do this, recalling the 
previous changes of variable, we further change the direction of the time axis in 
Eq. (|15[) and we deal with the following problem 

„ _ (1 + w 2 + w 12 ) [(2 - p)(l + w 2 ) - 2(p - 1) W] w' 
W (1 + w 2 ){l + w 2 + (p- l)w' 2 ) W ' [ ' 

under initial condition 

w(0) = 0, io'(0) = -a, a > 0. (43) 
We have the following result. 
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Proposition 5 Let p £ (1,2). There exists cyq > such that for any a < Qo, there 
exists a unique global solution to the problem (|42p - (|43p . and lim w(t) = 0. Further- 

more, the critical points ofw, that is, the solutions ofw'(t) = 0, consist of a sequence 
(in) such that w{t,2n) < and (t2n) ire local minima, w(t2n+i) > and (i2n+i) 
are local maxima. Moreover, we have iu(t2n— l) > wfon+l) for any positive integer n, 
and there exists no > such that \w{t n )\ is decreasing for n > no- In particular, the 
minimum is attained at t = to. 

Before proving Proposition [5] we transform (|42[) into an autonomous system 
whose associated phase plane will be analyzed in order to classify all the solutions. 
We divide this analysis into several steps. 

Step 1. Local analysis of the phase plane. We pass again to the autonomous 
system associated to (|42|) 

w = k, 

,_ (1 + w 2 + k 2 )[(2-p){l + w 2 )-2(p-l)wk]k ( 44 ) 
. k ~ (l + ™ 2 )(l + u; 2 + (p-l)fc 2 ) W ' 

Analyzing the phase plane, we find that (|44|) has a unique critical point in the 
plane, the origin, and the linearization near the origin has the matrix 

C(0 ' 0)= (- 1 -(2-P))' 

which has two complex eigenvalues with real parts (p — 2)/2 < 0. Thus, the origin is 
a stable node (attractor). From standard theory (see [211 Section 2.9]), we obtain 
that all the orbits of the system that lie in a neighborhood B(0,5) of the origin 
for some 5 > 0, will enter the origin in infinite time, meaning that 

lim w(t) = 0, lim w (i) = 0. 

t— too t—>oc 

Indeed, if w is defined on a maximal interval [0,T) with T < oo, then w can be 
extended to t = T in the C 1 class by simply letting w(T) = 0, and then continued 
afterwards, in contradiction with the maximality of the interval [0, T), whence 
T = oo. 

In order to establish the connections in the plane, we have to study the critical 
points of the system (|44[) at infinity. Since the system is not polynomial, we proceed 
as in the analysis in Section [3] and we find that the critical points at infinity are 
given by the polar angles 

I 7T 37T 3lT 77T / » — 1 \ / V — 1 \| , . 

fl£ < 7T, -jr-, —r, — , arccos = , arccos = > . (45) 

I 2 2 4 4 \ s /p*-2p + 2j' \^p 2 -2p + 2) J 

First of all, the local behavior near the nonhyperbolic critical points with polar 
angles 8 = ty/2 and 8 = 3ir/2 has already been studied in Section [3] We deduced 
that they are elliptic critical point. We denote by P3 the point with 8 = 3tt/2 and 
P3 the one with 8 = tv/2 in the sequel. 
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To study the local behavior near the other critical points, that are hyperbolic, 
we pass again to the Poincare sphere by letting 

' cos<P, k = — — — sin^. 



1-g 1-g 

After straightforward calculations, we obtain the following system 

' , _ e(i-e)[(i-g) 2 +e 2 ][2(p-i)e 2 cos<f sin*-(2-p)((i-e) 2 +e 2 cos 2 *)] sin 2 * 
8 ~ [(i-e) 2 +e 2 cos 2 *J[(i-e) 2 +-? 2 cos 2 *+(p-i)e 2 sin 2 *j ' 

< 

j:/ [(l-e) 2 + e 2 ] [2(p— l)e 2 cos* sin #-(2— p)((l-g) 2 + e 2 cos 2 * )] sin* cos* 

, — [(i-e) 2 +e 2 cos* 2 *J[(i-e) 2 +e 2 cos 2 *+(p-i)-j 2 sin 2 *j 



(46) 



In this setting, our critical points at infinity become Pi = (1, 7-7r/4), P2 = (1, arccos((p- 
1) / \/p 2 — 2p + 2)) and their symmetries P2 = (1, arccos(— (p — l)/- v /p 2 — 2p + 2)) 
and P* = (l,37r/4). Due to the symmetry of the system (|46p . we will analyze 
only the points Pi and P2. The linearization of (|46[) near the two points has the 
matrices 

M (l, l^j = I 2^ ° 2 



v 



p- 1 \ I 

AI 1, arccos — 



^P 2 -2p + 2 



1 



2 p--2p+2 
p P-1 



It follows that both points Pi and P2 are saddle points, of different orientation: 
there exists a unique orbit going out of Pi into the plane, while all the other orbits, 
that enter Pi, are contained in the curve g = 1, and there exists a unique orbit 
entering P2 and coming from the plane, the rest of orbits that go out from P2 are 
contained in the curve g = 1 and are not of interest to us. 

Step 2. No blow-up at the critical point P2. In this step we prove two prelim- 
inary facts that are needed before performing the global analysis of the previous 
phase plane. More precisely, we claim the following: 

(1) . There exists a unique positive, non-constant and decreasing solution to (|42|) 
corresponding to the orbit of the system (|46[) entering the critical point P2 and 
coming from the plane. 

(2) . The solution corresponding to the orbit entering P2 does this in infinite time. 
In order to prove assertion (1), we define 

g(w) : 



,'(t(w))- 

Following similar ideas and calculations to those in the proof of Proposition |H we 
deduce that g satisfies the equation 

g'(w) = i (g(w) 3 + (2-p)g(w) 2 + g(w)-F(w,g{w))} , (47) 

where 

2 g [-2(p -1X1 + w 2 )g 2 + (2 - p) 2 (l + w^g - 2(p - l)^ 2 ] 



F(w,g) 



(1 + w 2 ) [(1 + w 2 )g 2 + (p - l)w 2 } 
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Assuming that there are two different solutions w\, w 2 , corresponding to solutions 
.91, 32 of (147|) . we follow the same steps in the proof of Proposition [4] and by 
straightforward calculations, we obtain 

(.91 - 92)' {w) 



{91 - 92)(w) 



1_ 

w 



9iH + —^) + \92(w) + —^) + (gig 2 )(w) + l- 



(48) 



1 F(w, gi (w))-F{w,g 2 (w)) > g l FKfliH)-fKg 2 W) 
w gi(w) - g 2 (w) "I" w ffi(w)-32(w) 

for w sufficiently large. Recalling that we deal with orbits entering P%, we have 
lim g{w) = p — 1, whence the constant C satisfies 

C~^ + (p-1) 2 + 1-^^=p 2 . (49) 

On the other hand, by differentiating F(w, g(w)) with respect to its second variable 
g and passing to the limit as w — > 00, we easily get that F(w, g(w)) admits a 
Lipschitz constant for g for w large, which is exactly 2(p — 1). We skip the proof 
since it is similar to that of Step 4 in the proof of Proposition^ Since p 2 > 2(p — 1) 
for all p > 1, reasoning as there, for any e > 0, there exists wq = wo(s) sufficiently 
large such that 

gi{w)-g 2 (w) > ( w \° , (50) 
91 (too) - 32(^0) ~ \ w J 

where C = p 2 — 2p + 2 — 2e > for some e > small. Fixing such e > and wo, 
this implies 

lim (gi(w) - g 2 (w)) = 00, 

w—> 00 

in contradiction with the fact that lim gAw) = p — 1 for j = 1, 2. 

w— >oo 

In order to prove assertion (2), we introduce the following energy associated 
naturally to (|42)l : 

:= w 2 + w' 2 . (51) 

Thus, 

d _ 2(1 + w 2 + w' 2 )w' 2 [-(2 - p)(l + w 2 ) + 2(p - 1)W] 

dt ( l + ^2 )(1 + u; 2 + (p _ 1)w , 2) ■ ^ 

The idea is to bound from above the time derivative of E(w) by some function of 
E(w). Recalling the notation w' = k and the fact that w/k — > p — 1 for the orbit 
entering P 2 , we have 

d v( U\\ ot,±W< ^" ^ h(2^p)(l + ^ 2 ) + 2(p-l) W fc] 
Jt E(w(t)) = 2(1 + E(w(t))) {1 + w2){1 + w2 + {p _ 1)k2) 

= 2{1 + E{w(t)))H(w{t)). 

Estimating H(w) as w/k — > p — 1, we find that 

lim H(w) = — — - , 

W— >OD p — 1 
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hence there exists a sufficiently large constant C > such that 

j t E{w{t))<C{l + E{w{t))), t>to (53) 

for some to > 0. By direct integration in (|53[) . we deduce that the solution w is 
global. 

Step 3. Global analysis of the phase plane. In this step we end the analysis of 
the phase plane by establishing the connections between critical points. We restrict 
ourselves by symmetry to critical points in the hyperplane w < 0. 

As we already know, there is a unique orbit going out of Pi into the plane and 
a unique orbit entering Pj from the plane. We show that the orbit coming from 
Pi is attracted by the origin. Assume that this does not happen, thus this orbit 
may end up either at the elliptic point (l,37r/2) or at P2. In the latter case, since 
the phase plane has no self-intersections at finite points (by standard uniqueness 
arguments), by uniqueness of the orbits in Pi and P2, the orbits entering have 
to come from the elliptic point P3 = (l,37r/2), that will intersect the connection 
between Pi and P2 , contradiction. Finally, if the orbit from Pi goes to the elliptic 
point P3, the plane would remain incomplete, since the remaining two points are 
both attracting orbits from the plane. Hence, it remains the unique possibility 
that Pi is connected to the origin, and this is a global solution. Then, the plane 
is completed by the fact that the orbits from the elliptic point P3 go one of them 
to P2 and the rest either to the origin or back to the elliptic point P3, as shown 
in Figure[5] We will let then, ao := —w'(to), where w connects P3 and P2 and to 
is such that w(to) = and w(t) < for t > to- 

Thus, there are several global orbits. One of them is the one that connects 
Pi and O which (before the first intersection with the axe w = 0) corresponds to 
the unique minimizer. Coming back to h, we find that h(0) = 0, h(t) — > 7t/2 as 
t — > 00. We will study its oscillatory properties in the sequel. All the other global 
orbits come from the elliptic point P3 and reach it in a countable infinite number 
of times either forever or before being attracted by P2 or the origin. In any case, 
when coming back to h, all of them verify that lim h(r) = —00. These connections 

?'-s-0+ 

can be seen better in Figures [5] and [6] below. 



Remark 1 From the previous phase plane analysis, we deduce that 

(a) Let h* € C°°((0, 1]) be the solution to © with h*(l) = § , as given by Theorem 
03 We pass to logarithmic coordinates by letting f(t) := ft*(e*), t € (0,oo) and 
w{t) := cot/(i). Then, > w'(Q) > -a . 

(b) If w is a solution to (|42[) -(|43 p with a = ao, then w is a global solution and 
belongs to the orbit entering P2 in the phase plane. 

(c) If w is a solution to (|42p - (|4"3"l) with a > ao, then there is ti such that lim (w(t),w'(t)) = 

t—tti 

P 3 and w(t) < 0, w'(t) < for t G (0,ti). 

In order to characterize all global solutions to Q, we have, in a preliminary 
step, to exclude both that the equation becomes singular or that there is a blow-up 
in finite time. 
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Phase portrait for 1 <p<2 
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Fig. 5 Global phase portrait of 1441 . Numerical simulation for p = 1 + Note that all 
phase portraits are in fact symmetric with respect to w' . However, for clarity reasons, we plot 
only the orbits joining points in the lower half-plane plus one orbit in the upper half-plane. 



Phase portrait for 1 <p<2, zoomed at (0,0) 
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Fig. 6 Phase portrait of l]441 i near the origin. Numerical simulation for p = 1 + 

Lemma 4 Let h be a nonconstant solution of ([9|) in some interval (ro, ri) C (0, +oo). 
Then, 

(i) ( p -l)h'(r) 2 + s -^^>0forallre[r ,r 1 }. 
(ii) \h'(r)\ < C for all r € [r 0) n]. 

(Hi) h can be globally extended in (0, +oo) as a solution to ([9]). 

Proof Obviously, (iii) is implied by standard ODE theory by (i) and (ii). Moreover, 
both statements are equivalent stated in terms of nonconstant solutions in (£o,ti) 
of 

f"((p - l)/' 2 + sin 2 /) + (2 -p)/'(/' 2 + sin 2 /) = sin / cos /((3 -p)/' 2 + sin 2 /) (54) 
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via the smooth transformation /(f) = /i(e'): 

(a) f'(t) 2 + sin 2 (/(t)) > for all t G [t ,ti]. 

(b) \f'(t)\ < C for allie [to.ti]. 

Suppose that / is a solution of ()54[) in some finite time interval (to,ti). Then, 
performing the change of variables w = cot(/), we see that the only points where 

(a) or (b) may fail are to in case lim (w(t),w' (i)) = (+oo,— oo) and the points 

*->*+ 

T G (to, ii] such that lim (w{t),w'(t)) = (— oo,— oo), i.e. when the orbit of w in 

t-yT- 

the phase portrait approaches the elliptic critical point P3. 

Suppose now, by contradiction, that lim /'(t) 2 + sin 2 /(f) = (respectively 
lim f'(t) 2 + sin 2 f(t) = 0). Then, since at the elliptic point P3 we have 

lim ^imm = lim z™w 

t-vT- sm/(t) t-^T- ™(^j 

or, respectively, 

lim «»/(*)/;(*) = hnr Z# 
t->t+ am f{t) t ^t+ wit) 

we obtain that sin/(i) = o(/'(f)) in a left neighborhood of T (respectively in a 
right neighborhood of to). Then, substituting this in (|54[) . we obtain that 

f"(t) = t^.f'(t) + o(f'(t)), 

in a left neighborhood of T. By a simple integration, we find that f'(t) ~ K\e~ Kt 
as t -> T _ (resp. t — s> tj), for suitable constants K > 0, _Ki € K. This is a 
contradiction since /'(£) — >■ as i — >• T _ (resp. t — ¥ ij). Therefore, (a) is proved. 

For (6), it suffices to note that the right hand side of (|54p is negative when 
t — > T~ and the second member of the left hand side is positive when t — > T~ . 
Therefore, / must be concave in a left neighborhood of T, which directly yields 

(b) . In case lim \f'(i)\ = +00, then, as in case (a), one would obtain that 

/"(*) = ^/'(*) + *(/'(*)). 

in a right neighborhood of to- This leads to a contradiction if to is finite and thus, 
(b) is also proved in this case. 

Proof (Proposition 0) The proof is divided into two steps in order to ease the 
reading. We recall that we refer only to solutions w belonging to orbits that go to 
the origin O(0, 0) in the phase plane. 

Step 1. Oscillatory properties. In order to study the oscillatory properties of w, 
we need the following preliminary results 

• (A) Given a point s* such that w(s*) = and w (s*) ^ 0, there exists t* G 
(s*,oo) such that «/(t*) = and w(t*) < if u/(s*) < 0, respectively «/(t*) = 
and w(t*) > if w (s*) > 0. Moreover, t* is a local minimum (respectively local 
maximum) point for w. 



= —00, 



= +00, 
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• (B) Given a point t* such that w'(t*) = and w(t*) ^ 0, there exists 
s* G (i*, oo) such that w(s*) = and «/(s*) > if tu(t*) < 0, respectively w(s*) = 
and w'(s*) < if tu(t») > 0. 

Assume for the moment that claims (A) and (B) are true, and we follow with 
the proof of the Proposition. We apply (A) and (B) in an iterated form starting 
from s* = so = 0, where we know from (|43p that w(0) = 0, w (0) = —a < 0. We 
thus find an increasing sequence (t n ) such that w (tn) = 0, w(t2n) < 0, «)(i2n+i) > 
0, and, due to (A), (t2n) are local minima and (£271+1) are local maxima. Moreover, 
w is monotone on (£n,£ n +i) and there exists a point s n £ (£n>£n+l) such that 
u)(s n ) = 0, as (B) implies easily. The fact that £„ — > oo is easy: if t n — > T < oo, 
then by continuity of w and w' we get w(T) = w' '(T) = 0, thus w = by general 
theory, which is a contradiction. 

The fact that the local maxima are ordered, that is, iu(£2n+i) < w(t2n— 1)> 
follows directly from the ordering of the intersections of the orbit connecting Pi 
and the origin with the axis w = 0. In particular, the minimum is taken at £ = to 
as stated. 

Finally, for the last assertion, we recall the energy E(w) introduced in (|51[) . 
The sign of its derivative (|52[) is determined by the sign of the quantity 



G{w) := (2 - p)(l + w z ) - 2(p - 1) 



tow; 



From the phase plane analysis, we have that lim G(w(t)) = 2 — p > 0, hence 

t— 

there exists to > sufficiently large such that G(w(t)) > for t € (to,oo). Thus, 
E(w(t)) is decreasing on (to,oo). Let no be such that t n > to for n > no. Then 
E(t n ) > E(t n +i), which implies 

\w(t n )\ > \w(t n+ i)\, for any n > no, 

ending the proof. 

Step 2. Proof of statements (A) and (B). The proof of claim (A) is immediate 
from the phase plane. Indeed, if we have s* such that w (s*) < and w(s*) = 
and there is no t* £ (s*,oo) such that w (t*) = 0, it follows that to < forever 
and the orbit of w either goes directly to O(0,0) or it goes to P2 or P3. The latter 
is impossible since we deal only with orbits entering the origin, and the former is 
impossible since w' < implies w decreasing, contradiction. A similar argument is 
valid for s* such that w'(s*) > and w(s*) = 0. 

To prove assertion (B), we argue by contradiction. Assume we are at a point 
£* such that w(t*) < 0, w'(t*) = and there is no s* G (i*, 00) as in the statement 
of (B). Then, either the orbit of w goes directly to O(0,0) without cutting the 
axis w = 0, or it cuts again the axis w' = 0. The latter situation is eliminated in 
an obvious way, since after cutting again w = 0, w will become decreasing, thus 
it will self-intersect, which is not allowed in a phase plane. Suppose now that the 
orbit of w enters O(0,0) directly. We define 

m ■= % 

w(t) 

hence j(t*) = and j'(t*) = — 1, which is proved by noticing that w"(t*) = —w(t*). 
We then establish the equation satisfied by j, which is 

,_ .2 (1 + to 2 + w 2 j 2 ) [(2 — p)(l + w 2 ) — 2(p — l)w 2 j] j 

3 3 (l + w 2 )(l + w 2 + (p-l)w 2 j 2 ) • 1 ; 
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Thus, in a sufficiently small neighborhood of the origin we have 



•/ -i -2 

J < "I - J , 



It follows that j blows up to — oo at some point s* > £*. Hence w(s*) = and, 
by the assumption, «/(s*) = 0. Coming back to /, this implies that f(t) = ty/2 
for t > s*. On the other hand, we know that / is nonconstant, thus, by part (i) 
in Lemma 21 (|54|) is nonsingular and we can apply standard uniqueness theory to 
get / = 7r/2, contradiction. This finishes the proof. 

We are now ready to prove Theorem [2] 

Proof (Theorem^) 

(a) and (b) are a direct consequence of Remark [1] and Proposition [5] 

To prove (c), for £ < ^ we scale r so that h(f) = h(ar) satisfies /i(l) = £ with 
a > 0. By the scale invariance of (O, h is also a solution, and h is increasing and 
positive in (0,1). Hence, by Proposition [3] it coincides with the minimizer hi of 
Gp. 

To prove (d), we have to characterize finite global non-constant solutions h to 
Q. Take one of them: we may assume without loss of generality that h(r) £ (0,7r) 
for r £ (ri,r2) C (0,+oo). Then, we can pass again to (|42p with to— variable in 
some time interval (ti,ta). By the analysis of the phase plane, since h is finite and 
global, then w{t) = cot/i(e*) is such that (w(t),(w)'(t)) joins Pi and the origin. 
Since this orbit is unique, then by the invariance of Q, h is of the form in (d). 

Finally, again by the analysis of the phase plane, there are infinitely many orbits 
connecting the elliptic point P3 and the origin and a unique orbit connecting P3 
and P2 . When translating this to the original variable h, by Lemma 3] all these 
solutions can be globally extended and a direct computation shows that h(0) = —00 
or h(0) = +00. In the former case, again by the invariances of (J9j> , h is globally 
increasing and has either +00 or kir as a limit for some feeZ (cases P3 and P2) 
or < ~ 2k + 1 ^ for some k G Z (in case the orbit joins P3 and the origin). The latter 
case is totally symmetric. 

5.2 The case p > 2 

In this subsection we prove Theorem [3] We start our study from Eq. (|42p and 
the associated autonomous system (|44[) . Analyzing the phase plane, we find that 
the system has the same critical points as in Subsection 15.11 the origin in the 
plane and the six different critical points at infinity, corresponding to polar an- 
gles given in J45[) . more precisely, on the Poincare sphere, Pi = (l,77r/4), P2 = 
(1, arccos((p— l)/\/p 2 — 2p+ 2)) and their symmetric points P2 = (l,arccos(— (p— 
1) /\/p 2 — 2p + 2)) and Pf = (l,37r/4), together with the elliptic critical points 
P3 = (l,37r/2) and P3 = (l,7r/2). By symmetry, we analyze only the lower half- 
plane. 

The linearization of the system (|44[) near the origin has the same matrix 
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which has two complex eigenvalues with real parts (p— 2) /2 > 0. Thus, the origin is 
now an unstable node. This is the only difference for the local analysis with respect 
to the case 1 < p < 2, studied in Subsection l5.ll We recall that the critical points 
at infinity, whose analysis is identical, are: P± and P2 saddle points of opposite 
orientation, that means, from Pi there is a unique orbit going out into the plane 
and there is a unique orbit entering P2 from the plane. With all these, we are able 
to perform the global analysis of the phase plane. 

Global analysis of the phase plane. There is a unique orbit going out from Pi 
into the plane. This orbit could enter the elliptic point P3 or P2. Assume by con- 
tradiction that there exists an orbit connecting Pi and P%. Then, the orbits going 
out of the unstable node 0(0,0) cannot go to the elliptic point P3, as they would 
intersect the orbit connecting Pi and P2, nor to P2, as P2 admits a unique orbit 
coming from the interior of the plane, that would be in our assumption the one 
coming from Pi. But this means that 0(0,0) will remain isolated, contradiction. 
Thus, the orbit coming from Pi enters the elliptic point P3. The orbits starting 
at O(0,0) will thus enter P3 too, except one of them which will go to P2. We 
can thus classify all the orbits of the system (|44[) into the following four types of 
connections: 

• A unique orbit coming from Pi and connecting P3. Going back to h, it gives rise 
to an increasing solution with h(0+) = 0, and h(ro) = -it at some ro > 0. As we 
shall see later, this orbit will be extended beyond ro into a global solution with 
h(r) — > 00 as r — > 00. 

• An infinite number of orbits starting at O(0, 0) and entering the elliptic point P3. 
Since, due to periodicity, we restrict our analysis to h(0+) € [0, 7r), going back to 
h, they represent solutions with h(0) = n/2. These solutions may oscillate around 
7r/2 in some finite interval (0,ro) and then increase to h(ri) = n for some ri > ro 
(or decrease to h(ri) = for some r\ > ro). 

• A unique orbit coming from O(0, 0) and entering P2. Going back to h, this orbit 
contains solutions with the following properties: h(0) = tt/2, h oscillates around 
7r/2 in some finite interval (0,ro) for some ro > and then it increases such that 
h(r) — > ty as r — > 00 (and by symmetry, there are also solutions that decrease with 
h(r) — > as r — > 00, the ones corresponding to the orbit connecting O(0, 0) to P%, 
the symmetric of P2 in the plane) . 

• Elliptic orbits going forever around the elliptic point P3. These orbits will contain 
solutions with h(0+) = — 00 and h(r) — > 00 as r — > 00. 

The connections described above can be seen better in Figures and [8] below. 

In order to characterize all global solutions to ([9]) for p > 2, we have to exclude 
again both that the equation becomes singular or that there is a blow-up in finite 
time. 

Lemma 5 Let h be a nonconstant solution of ([9|) in some interval (ro, ri) C (0, +00). 
Then, 

(1) (p-l)h'(r) 2 + s -^^>Oforallre[r ,n}. 
(it) \h'(r)\ < C for all r G [r ,n]. 

(Hi) h can be globally extended in (0,+oo) as a solution to ([9|). 
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Phase portrait for p>2 
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Fig. 7 Global phase portrait of H44II for p > 2. Numerical simulation for p = 3. 




Fig. 8 Phase portrait of 1441 near the origin for p > 2. Numerical simulation for p = 3. 



Proof We proceed as in the proof of Lemma[U First of all, (iii) is implied by (i) and 
(ii). Next, we pass to logarithmic coordinates by f(t) = h(e ) and to w = cot/, 
as in Lemma [3] (ii) is proved as (i) in Lemma [3] with the only difference that 
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supposing that lim f (t) = +00 at the elliptic point, we arrive at 
t->T- 

/"(*) = fff/'W + "(/'(*))- 

which is again a contradiction if T is finite. 

For (i), we note that it is equivalent to prove that /' does not vanish and, as 
in Lemma HI it suffices to check it at the elliptic points for the phase portrait of 
w. We perform now the following change of variables: g = Then, by (|42p . we 
have 

,_ w 6 g 2 ((p-2)wg + p-2-(p-l)g 2 ) 
9 w{l + w 2 )(l + w 2 +g 2 w 4 ) 

w i ((p-2)g(l + g 2 ) + (p~5)g 2 -l) 
w{l + w 2 )(l + w 2 + g 2 w 4 ) 

2w 2 ((p - 2)wg - 1 - g 2 ) + ((p - 2)wg - 1) 
w(l+w 2 )(l + w 2 +g 2 w 4 ) 

Noting now that w(t) — > —00 as t -> T~, g < in a left neighborhood of T, we 
see that if git) — > CP as t — > T~ , then, since w(t)g(t) = ~ ► +°° as f — T, we 
get that all terms in the right hand side of the above equation for g are negative. 
Therefore, g(t) -fr _ . This yields the conclusion, since 

Combining Lemma[S]and the global analysis made before, we rigorously deduce 
that the orbits going to the elliptic point P3 contain indeed global solutions of ([9]) , 
as we already anticipated when doing the analysis. 

The proof of Theorem [3] is analogous to that one of Theorem [2] and we skip it. 

Remark. The special case p = 2. The case p = 2 has been widely studied with 
the use of direct methods (see e.g. the referenced works [21H]). In this case, there 
exists a unique minimizer for the energy and it has (in our notations) the explicit 
formula 

h(r) = 2 arctan r. 

This particular case can be easily understood with our technique of phase plane 
analysis, as it matches perfectly between the pictures for 1 < p < 2 and for p > 2 
respectively. We complete thus the picture of the evolution of the solutions with 
respect to the variation of p. We start again from the system (|44p and study the 
phase plane. The critical points are the same O, Pi, P2, P3, P*, P£ , P3 as above, 
and the local analysis around the critical points at infinity remains identical. In 
change, it can be shown that O(0, 0) becomes a center. Thus, passing to the global 
analysis of the phase plane associated to the system (|44[) . and restricting our 
analysis to the lower half-plane, we notice that: 

• Exactly in this case, the unique orbit going out from P\ into the plane, enters 
P-2- This orbit contains the solution h(r) = 2 arctan r. 

• All the rest of orbits are either elliptic orbits around P3, or orbits that oscillate 
forever around the center O(0, 0). The latter orbits contain solutions which oscillate 
forever around it/2. 
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The connections described above can be seen better in Figures l9l and [TOl below. 
Indeed, the connection Pi — P2 is the special and explicit one in this case (and, as 
we have seen, it only exists in this case) . 
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Phase portrait for p=2 



-101 





Oscillatory orbit 


Orbit connecting 

P1 and P2 .' 

.> /! 

1 • 
1 ■' 


H X 

■ 1 

'- = *. 

1 ' \ 

s*6 Elliptic orbits 

■ 1 

■ - : 

1 



00 



100 



Fig. 9 Global phase portrait of i ]44tl . Numerical simulation for p = 2. 
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Phase portrait for p=2, zoomed at (0,0) 
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Fig. 10 Phase portrait of 1144 II zoomed near the origin for p = 2. 
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